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Abstract: A competition model on N^ between three clusters and governed by di- 
rected last passage percolation is considered. We prove that coexistence, i.e. the 
n~! i three clusters are simultaneously unbounded, occurs with probability 6-8 log 2. 

When this happens, we also prove that the central cluster almost surely has a pos- 
itive density on N^. Our results rely on three couplings, allowing to link the com- 
petition interfaces (which represent the borderlines between the clusters) to some 
particles in the multi-TASEP, and on recent results about collision in the multi- 
TASEP. 

> 
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1 Introduction 



The directed last passage percolation (LP?) model has been much studied recently. 
In dimension 2, it is closely related to some queueing networks, to random matrix 
C^ ' theory and to some combinatorial problems such as the longest increasing subse- 

quence of a random permutation. See Martin [9 | for a quite complete survey. 

Throughout this paper, N denotes the nonnegative integer set. We consider 
i.i.d. random variables a>(z), z e N^, exponentially distributed with parameter 1. 
Let P be the Borel probability measure induced by these variables on the product 
space Q. = [0, oo)^ . The last passage time to z is defined by 

G{z) = max V w(z') 

z'ey 
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where the above maximum is taken over all directed paths from the origin to z (see 
Section[2]for precise definitions). The maximum G{z) is a.s. reached by only one 
path, called the geodesic to z- As a directed path, this geodesic goes through one 
and only one of the three sites (0, 2), (1,1) and (2, 0), called sources. Let the cluster 
C{s) be the set of sites z e N^ whose geodesic goes by the source s. Hence each 



configuration oj 
Figure [1] 



e Q yields a random partition of {{x,y) e W : x + y > 2\, see 





Figure 1: Two simulations of the clusters C(0, 2), C(l, 1) and C(2,0) which have 
been colored with respectively dark blue, light blue and red. To the left, C(l, 1) 
seems to be unbounded (there might be coexistence) whereas, to the right, it is 
bounded. Note that such a simulation of bounded but large cluster C(l, 1) is very 



rare. 



We focus on the competition (in space) between the three clusters C(0, 2), C(l, 1) 
and C(2, 0). The directed character of the model implies the first and the third ones 
are unbounded. But this is not necessary the case for the second one; we will talk 
about coexistence when the cluster C(l, 1) is unbounded. 

Our main result (Theorem [T]) states that coexistence occurs with probability 
6-8 log 2, which is close to 0.4548. As far as we know, there is no other model 
where such a coexistence probability is exactly computed. For instance, in the 
(undirected) first passage percolation model, the competition between two clusters 
growing in the same space leads to two situations: either one cluster surrounds 
the other one, stops it and then infects all the other sites of 1? or the two clusters 
grow mutually unboundedly, which is also called coexistence. And in the case 
of independent exponential weights, Haggstrom and Pemantle 171 have proved 
that coexistence occurs with positive probability. Caret and Marchand @ have 
since generalized this result to ergodic stationary passage times and to random 
environment. 
Our second result (Theorem |2]) completes the first one. When the cluster C(l, 1) is 



unbounded then it almost surely has a positive density in the following sense: 



limsup — Card(C(l, 1) n [0,nf) > , 



The proofs of Theorems [T] and [2] are mainly based on three couplings: see Tho- 
RissoN [ 13] for a complete reference on couplings. The first one is due to Rost. In 
IfTTl . he builds a totally asymmetric simple exclusion process (TASEP) from the 
LPP model, using the last passage times G(z), z e N^, as jump times. A back- 
ground on exlusion processes can be found in the book IS] (Part III) of Liggett. 
The borderlines between the clusters C(l, 1) and C(0, 2) and between C(l, 1) and 
C(2, 0) are modeled by two infinite directed paths, called the competition inter- 
faces. Ferrari and Pimentel [5J, thus Ferrari, Martin et al [4] have studied their 
asymptotic behaviors. These competition interfaces play an important role here 
since the cluster C(l, 1) is bounded whenever they collide. The Rost's coupling 
allows to link these competition interfaces to two tagged pairs [oo 1] in the TASEP, 
where labels oo and 1 respectively represent holes and particles. In particular, the 
coexistence phenomenon is equivalent to the fact that these two tagged pairs never 
collide (Lemma 13. 

The second coupling allows to turn the two tagged pairs into two second class par- 
ticles whose labels are denoted by 2 and 3 (Lemma|3]l. A second class particle is an 
extra particle which interacts with particles like a hole and interacts with holes like 
a particle. Its trajectory has been studied by Mountford and Guiol flO|. See also 
Seppalainen fn\. The idea to represent a second class particle as a hole-particle 
pair [cxi 1] is due to Ferrari and Pimentel 151 . 

Ferrari, GoNfALVES et al Q have studied the collision phenomenon of two sec- 
ond class particles. Thanks to the two previously announced couplings, they de- 
duced (Theorem 4.1) that coexistence occurs in the LPP model with probability 
1/3. However, they assume for that some constraining initial conditions, namely 
aj(0, 0) - a»(l,0) - a>(0, 1) = 0. We will explain why their coexistence result is a 
partial version of Theorem [T] 

Finally, the third coupling, usually called basic coupling (|'8l, p. 215), allows to 
consider the two second class particles (i.e. the 2 and 3 particles) in a more general 
exclusion process, the multi-TASEP . Recently, Amir et al HI have proved many 
results about this process. Some of them are expressed in terms of second class 
particles (Proposition |4] and Lemma|5]l, thanks to that third coupling. 
To sum up, these three couplings state a strong link between the multi-TASEP and 
the LPP model, leading to Theorems [T] and |2] 

The paper is organized as follows. Section |2] contains the definition of the LPP 
model and statements of main results with some comments. Section lXTl introduces 
the TASEP. The tagged pairs [oo 1] are identified in Section IX2l Sections [3.31 and 
13.41 are respectively devoted to the second and the third coupling on which proofs 
of Theorems [T] and |2] are based. The first coupling is described in Section 14.11 
Competition interfaces are defined in Section 143] and linked to tagged pairs in the 



TASEP in Section l44l Finally, Theorems [T] and [2] are proved in Section |5] 

2 Coexistence results 

Recall that P denotes the law on O = [0, oo)^ of the family {coiz), z e N^) of i.i.d. 
random variables exponentially distributed with parameter 1. 
A directed path y from (0, 0) to z is a finite sequence of sites {zq,z\,- ■ ■ , Zk) with 
zo = (0, 0), Zk = z and z,+i - Zt = (1, 0) or (0, 1), for < / < A: - 1. The quantity 
Yjz'ey ^(z') represents the time to reach z via y. The set of all directed paths from 
(0, 0) to z is denoted by r(z). The last passage time to z is defined by 



G{z) - max V w(z') 



z'ay 

Since each path of r(z) goes by either z - (1 , 0) or z - (0, 1), the function G satisfies 
the recurrence relation 

G{z) = coiz) + max|G(z - (1,0)), G{z - (0, 1))) (1) 

(with boundary conditions G(z) = for z = (x, -1) or (-1, x) with x e N). A site 
z is said infected at time t if G(z) < t. Relation ([T]) can be interpreted as follows: 
once both sites z - (1, 0) and z - (0, 1) are infected, z gets infected at rate 1. 

Recall that the cluster C{s) is the set of sites z € N^ whose geodesic goes by 
the source s. Let us point out the directed character of the LPP model forces the 
clusters C(2,0) and C(0, 2) to be unbounded. Indeed, if the site z = {x,y) belongs 
to C(2, 0), so do all the sites on its right. Similarly, if the site z = {x,y) belongs to 
C(0, 2) so do all the sites above. Actually, only C(l, 1) can be bounded. Indeed, 
whenever 

min{w(l, 0) + cu(2, 0), oj(0, 1) + oj(0, 2)) > w(l, 1) + max{w(l, 0), w(0, 1)} , (2) 

the last passage times G(2, 0) and G(0, 2) are both larger than G(l, 1). In this case, 
sites (2, 1) and (1,2) respectively belong to C(2, 0) and C(0, 2), hence the cluster 
C(l, 1) is reduced to its source. See also the right hand side of Figure [T] for the 
simulation of a larger (but bounded) cluster C(l, 1). 
For any positive integer n, let 

a(?i)-Card(C(l,l)n{(x,3;)eN^ : x + y = n}) . 

We will say there is coexistence when the cluster C(l, 1) is unbounded, i.e. a{n) > 
for all n >2. When this holds, each cluster C(s) contains sites whose last passage 
time is as large as wanted; the three clusters C(0, 2), C(l, 1) and C(2, 0) coexist. 

Theorem 1. Coexistence probability is 6 - S log 2: 

P(V« > 2, a{n) > 0) - 6 - 8 log 2 . 



It is already known that coexistence probability differs from and 1 . Indeed, 
it is clear that coexistence cannot hold a.s. since the event defined in Q occurs 
with positive probability. Moreover, in a previous work [2|, we have shown in 
particular that coexistence occurs with positive probability if and only if there exists 
at least one infinite geodesic (different from the horizontal and the vertical axes) 
with positive probability; this last condition being proved in [5 ], Proposition 7. 

Let us compare our result to Theorem 4.1 of |3|. In that paper, Ferrari, 
GoN^ALVES et al prove that coexistence occurs with probability 1/3, but they con- 
sider the LPP model under the initial condition 

w(0,0) = a>(l,0) -w(0, 1) = 0. (3) 

Since the origin (0, 0) belongs to each geodesic, its weight does not affect the coex- 
istence probability. However, the cluster C(l, 1) benefits from max{c(;(l, 0), a»(0, 1)) 
whereas the clusters C(2,0) and C(0, 2) only use a»(l,0) and a>(0, 1) respectively. 
Assuming a>(l, 0) = a>(0, 1) = amounts to remove this benefit. More precisely, let 
g : Q ^ Q defined by g(w)(0,0) - ^(w)(l,0) - g(w)(0, 1) = and g(w)(z) - aj{z) 
otherwise. It then follows 

C{\,l){g{aj))cC{\,l){oj). 

Theorem 4.1 of fT\ says C(l, l)ig{a))) is unbounded with probability 1/3. This 
suggests that coexistence probability in the LPP model (without initial conditions) 
is greater than 1/3. Actually, this remark has motivated the present work. 

Our second result concerns the density of the cluster C(l, 1) in the quadrant 
N^. Let us first remark that if the density of the cluster C(l, 1) is positive, i.e. 

lim sup —Card (C(l, 1) n [0, nf) > , (4) 

then C(l, 1) is unbounded. Moreover, (01) holds if and only if 

1- '^<^") ^ A 
lim sup > 0. 

This stems from the fact that a{n + I) belongs to {a{n) - 1, a{n), a{n) + I], for any 
n. Hence, the inequality a{n) > 6n for some 5 > and integer n, implies that the 
cluster C(l, 1) n [0, n]^ contains a square with diagonal of length \6n\ (where \_x\ 
denotes the integer part of x). 



Theorem 2. Coexistence almost surely implies positive density for C(l, 1)." 

Sn > 2, a{n) > 1 = 1 . 



„ ,,. C({n) 

P I hm sup --- > 

rt— >co n 



Moreover, 



a(n) 
P(limsup— -- < 1| = 1 . 



3 TASEP and related processes 

3.1 Some definitions 

In the sequel, TASEP stands for totally asymmetric simple exclusion process. It 
is a Markov process whose dynamics can be easily described: at rate 1 (i.e. after 
an exponential time with parameter 1), particles (integer or oo) at sites x and x+ \ 
attempt to exchange their positions. The exchange occurs if the value at site x is less 
than the value at site a + 1, otherwise nothing happens (total asymmetry property). 
There is at most one particle per site (exclusion condition). The oo particle has thus 
a role of hole. Here is a precise definition: 

2^ 

Definition 1. Set Z = Z U {oo}. Let S be a subset ofZ . Consider the linear 
operator X on cylinder functions f on S defined by 

xeZ 

where rj^'^'^^ is obtained from rj = {rjy,y e Z} by exchanging values at x and x + I: 



nr' 



rjy ifyi{x,x+l}, 
rjx+i ify = X, 
Vx ify = x+ 1. 



A Markov process on R+ with configuration (or state) space S and with generator 
X is called 

(a) TASEP if the configuration space is S = {1, oo)^, 

(b) k-type TASEP if the configuration space is S ={1,2,...,^, oo}^, 

(c) multi-TASEP if the configuration space is S = Z^. 

Let us add that the order relation < on Z is extended to Z as follows: r < oo if 
and only if / belongs to Z. 

Besides, it will be convenient to locate some particles of interest in a configuration. 
Let ?7 be a configuration in 5 c Z containing exactly one k particle (k e Z). The 
position of this k particle in t] is denoted by 

k[Ti]. (6) 

For a further use, it is convenient to introduce the following particular configura- 
tions, described in Figures [2l [3] and |4] For any integer m, 

• Let 7]'" e { 1, oo)^ defined by 

™^|l ifxe{...,-3,-2}U{0)U{m + 2), 
oo otherwise. 



Let T?^^)'"' e {1, 2, 3, oo)^^ defined by 

'l if xe {...,-3,-2,-1}, 



(3),m 
Wx = 



2 if X = 0, 

3 if X = m + 1, 
oo otlierwise. 



(8) 



• Let 77(") € Z^ defined by 

rlr^ = x (xeZ). (9) 

3.2 Tagged pairs in the TASEP 

We want to follow the evolution of two pairs of particles over time in the TASEP 
with initial configuration 77™ defined in (|7]). A pair consists of a couple (00, 1) 
tagged with brackets. In the configuration if, there are exactly two pairs [00 1], the 
left one is called - pair and the right one + pair (see Figure |2ll. 



4)] 



§ 



(X3^ (ex,) (1 



pair ' + pair 



Figure 2: Configuration 77'" with the two tagged pairs [00 1]. They are separated by 
m "holes" 00. On the axis Z, the origin is marked with a vertical arrow. 

Let us describe the evolution rule of the two pairs. Let e € {-, +). When a 1 
particle jumps (from the left and at rate 1) over the hole 00 of the e pair, this one 
moves one unit to the left: 

rv 

l[ool] becomes [ooljl. (10) 

When the 1 particle of the e pair jumps to the right (over a hole 00 and at rate 1) 
then the e pair moves one unit to the right: 

[ool]cx3 becomes oo[oo 1] . (11) 

Definition 2. For e e {-,+), let us denote by H^{t) the hole's position of the e pair, 
at time t, in the TASEP with initial configuration r/'": 

■ • • [ 00 1] • • . [ 00 1] • ■ • 

H-{t) H+(t) 

The collision time is defined as 

Teoi = infU > : H-(,t) = H\t)} , (12) 

with the convention inf % = oa. 



The two tagged pairs merge together at the jump time Tcoi and there remains 
only one tagged pair so that H~{t) - H^{t) =: H{t) for all t > T^oi'- 

[ CX3 1][ oo 1] becomes oo[ oo 1]1 . 

H-(t) HHt) Hit) 

Finally, let us point out that the process {H^it), t > 0), for e e {-, +), is not marko- 
vian but {{^{t),H~{t),H-'{t)),t > 0) is, where ^ = {^{t),t > 0} is a TASEP with 
initial configuration rf. The reader can refer to tSj to get more details on tagged 
particle processes. 



3.3 From tagged pairs to three-type TASEP 

Recall that in a /c-type TASEP (or in a multi-TASEP), a / particle can pass a j 
particle if and only if / < j. But the above evolution rule shows that each tagged 
pair behaves like any single / particle with respect to a oo particle -see (fTTT) - and 
also with respect to a 1 particle -see (ITOl)- provided / is more than 1 and finite. If 
we turn the - pair into a 2 particle and the + pair into a 3 particle (for instance), 
we obtain a three-type TASEP. More precisely, consider transformations 



mJ; 



(^?'\e2 



:{1,~} 



{1,2,3, ex,) 



defined by 














(a) For jc + 2 < J, 














Vz-l 


ifz 


< X, 
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ifz 


- x-l- 1, 






^r(v) = ■ 


Tlz 


iix + 2<z<y- 


-1 
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ifz 


- y, 








,'7z+i 


ifz 


>3'+ 1. 




(b) For X = y, 














Vz 


-1 


ifz<x-l. 






^?'('7) - 


3 

2 




ifz-x, 
ifz-x+1. 










,'7z+i 


if z > x + 2. 





For example,^"^''"''"^ transforms 77™ (Figure|2]l into 77^^^-'" (Figure[3ll. 

In what follows, we focus on the evolution of the two particles 2 and 3 over 
time in the three-type TASEP until the collision time Tcoi. The applications ^F^-^' 
provide the following coupling: 

Lemma 3. Let ^ - {^(t), t > 0] be a TASEP with initial configuration rj'" and 
collision time T^oi as defined in (1121 ). Then, the process 



f ■- ["¥"' 



{t),H*(t) 



(^(0), < f < Teol 



® (2 



(oo ) (kT) 



Figure 3: Configuration rp'^^'". The two particles 2 and 3 are separated by m holes 



is a three-type TASEP on [0, Tcoi] with initial configuration t;^^)''". In particular, 
with the notation ^, it follows 

(i) For t < Tcoi, 2[f (?)] = H-{t) + 1 and 3[f (?)] - H^{t), 

(ii) For t = Tcoi, 2[f (?)] = //""(O + 1 = //"(O + 1 and 3[f (?)] = //''(O - H-{t), 

It is crucial to remark this coupling holds until time Tcoi (Tcoi included thanks 
to the part (b) in the definition of ^'^■^ ). 
Particles 2 and 3 in the three-type TASEP ^' can be seen as second class particles. 

3.4 From three-type TASEP to multi-TASEP 

The goal of this section is to couple a three-type TASEP with initial configuration 
^(3),m ^j^j ^ multi-TASEP with initial configuration 77*^°°^ using the basic coupling 
(see lH). To do so, let us consider a family {Nxit), t > 0, x e Z] of independent 
Poisson processes with parameter 1 . At each event time Nxit) and for the two pro- 
cesses, the particles located respectively at site x and x+l exchange their positions 
if permitted by the order <, nothing changes otherwise. Hence, the two processes 
evolve simultaneously on the same probability space. See Figure |4l 



© 




(m+l) 673+2) 



Figure 4: The configurations 77^^^''" and 77^°°^ are the starting points of the three-type 
TASEP and the multi-TASEP under the basic coupling. They are aligned so that 
2 and 3 particles in the three-type TASEP respectively correpond to and m+l 
particles in the multi-TASEP (at time t = 0). 

First, let us remark some occurring jumps for the multi-TASEP, say between a 
/ particle and a j particle (with / < j), are not authorized for the three-type TASEP. 
This happens when the corresponding particles in the three-type TASEP are the 

9 



same or when / € { 1 , . . . , m) and the corresponding particle to j in the three-type 
TASEP is 3. Then, we deduce that up to the time where the 2 particle passes the 3 
one in the three-type TASEP, 

• the 2 particle in the three-type TASEP corresponds to the particle in the 
multi-TASEP; 

• the 3 particle in the three-type TASEP corresponds to the further right parti- 
cle among particles 1 , . . . , m -i- 1 in the multi-TASEP. 

Hence, the time where the 2 particle and the 3 particle exchange their positions in 
the three-type TASEP is also the time where the particle has just overtaken all 
the particles 1, . . . , m -i- 1 in the multi-TASEP. Theorem 7.1 of Amir et al. m states 
this last event occurs with probability 2/(m -i- 3). Now, the basic coupling allows to 
transfer this result to the three-type TASEP: 

Proposition 4. Let P'„ be the probability measure of a three-type TASEP £,' with 
initial configuration 7/*^^^'". V</ith notation Q, it follows 

p;(3?>o,2[r(o]>3[r(o]) 



m + 3 



Note that, before results of |[T1, this result had been conjectured (and proved in 
the case m e {0, 1)) by Ferrari et al in |[3l. 

Let us respectively denote by ^' and ^" a three-type TASEP and a multi-TASEP 
with initial configurations 7^(3),m ^^^^ ^(oo) Uj^^jj jj^g gj^^j ^f jj^j^ sggtion, we assume 

that Vf, 2[^'(f)] < 3[^'(f)]- The basic coupling described above implies, at any 
time t, the 2 particle in the three-type TASEP corresponds to the particle in the 
multi-TASEP, i.e. 

V? > 0, 2[^'{t)] = 0[r'(0] , 

and the 3 particle in the three-type TASEP eventually corresponds to one of the 
particles 1, . . . , m -i- 1 in the multi-TASEP, i.e. 

3ke{\,...,m + \},3tk,\ft>tk, 3[f (?)] - k[^"{t)] . 

The fundamental result (Corollary 1 .2) on which fP| is based is that in the multi- 
TASEP with initial configuration r]'^°°\ each particle chooses a speed. Precisely, for 
every k eZ, 

lim = Uk a.s. , 

where {Uic,k e Z) is a family of random variables, each uniformly distributed on 
[-1; 1], called the TASEP speed process. So, on the event {Vt,2[^'{t)] < 3[f (f)]), 
the ratios 2[^'{t)]/t and 3[^'{t)]/t converge respectively to Uq and Ut, for a given 
k. To sum up, the event 

r 3[f (Q] - 2[f (Q] H V. or^VAi ^ ^r^VAi 

lim and V?, 2[^ (?)] < 3[^ (f)] 

10 



is a.s. included in 

m+l 



\J [Uo = Uk and V?,0[r'(0] < k[^"m . (13) 

k=i 

Finally, Lemma 9.9 of H) states, in the multi-TASEP with initial configuration 
T]'-"'\ every two particles with the same speed swap eventually. So the event (IT3l) 
has zero probability. 

Lemma 5. Let P'„ be the probability measure of a three-type TASEP ^' with initial 
configuration 7/*^^^"'. Then, 

P'Jlim 3[^'(0] - 2[r(0] ^ Q ^^^ Vf,2[^'(0] < 3[^Xt)]] = . 

\f^oo t I 

4 LPP model and tagged TASEP 

The goal of this section is to state a coupling between the LPP model and the 
TASEP. This coupling allows to link the competition interfaces (defined in Section 
14.31) to some pairs of particles (identified in Section l42l) . 

4.1 Rost's coupling 

In ifTTTl . RosT gives an explicit construction of the TASEP from the LPP model, 
using the last passage times G{z), z e N^, as jump times. Let us describe this 
construction. 

Let us start with the configuration 77'''^' € {1, 00}^ which is made up of 1 particles 
on nonpositive integers and 00 particles on positive ones. The Rost's idea consists 
in labelling 1 particles from the right to the left hy Pq, Py, P2- ■ ■ and 00 particles 
from the left to the right by Hq, Hi, H2. . . as in Figure |5] and in following them 
over time. Letters P and H refer to particle and hole. 




Figure 5: Here is the configuration r]'' with labelled particles. On the axis Z, the 
origin is marked with the vertical arrow. 

The evolution rule is 

Pj and Hi exchange their positions at time G{i, j). (14) 

At time G(0, 0) - a»(0, 0) the first exchange takes place between Pq and Hq. The 
second one will concern Pq and Hi if w(l, 0) < a;(0, 1) and Pi and Hq otherwise. 



11 



More generally, at time max{G(/ - 1, j), G{i, j - 1)), the exchanges between Pj and 
Hi-i, and between Py_i and //,• have already taken place. Labels of 1 particles 
and those of oo particles remaining sorted over time then Pj is then the left nearest 
neighbor of //,. From that moment, they exchange their positions after the time 
a){i, j) (i.e. at rate 1) thanks to the recurrence relation ([T|l: 

co{i, j) = Gii, j) - max{G(/ - 1, j), G{i, j - 1)) . 

It then suffices to disregard labels Pj and //, to get back the TASEP. Precisely, let 
us denote by Pjit) and Hi{t) the positions of particles Pj and //,• at time t. At the 
beginning, P/0) = -j and Hi(0) = i + \. Now, set for ? > and x eZ 

II if there exists /' such that P At) - x, 
oo otherwise, 

and let ^{t) be the configuration i^xit))xez- Then: 

Lemma 6. The process ^ = {^(t), t > 0} is the TASEP with initial configuration 

^ext 

Let us end this section with describing an explicit way to obtain the configura- 
tion ^(t) from the infected region at time t, i.e. the set {z e N^ : G(z) < t]: 

1. In the dual lattice (-5,-5) + N^, draw the border of the infected region 
at time t and extend it on each side by two half-line, as in Figure |2l The 
obtained broken line consists of horizontal and vertical unit segments; it rep- 
resents the axis Z on which ^(t) is defined. 

2. Mark the last (from north to east) unit segment of the broken line before the 
diagonal y = x;it represents the origin of Z. 

3. Replace each vertical (resp. horizontal) unit segment of the broken line with 
a 1 (resp. 00) particle. 

For instance, the configuration of Figure [6] is obtained thanks to the previous algo- 
rithm from the infected region given by Figure |7J 

4.2 Initial conditions in the LPP model 

Consider the integer valued random variable N defined by 

max{m > 1 : a>{\,0) -1- ... -1- a>{m, 0) < a»(0, 1)) if exists. 



otherwise. 

We first remark that {A'^ > 1} = {w(l,0) < co(0, 1)} occurs with probability 1/2, by 
symmetry. 



Lemma 7. Conditionally to {N > 1 ), the random variable N is distributed accord- 

l 

2- 



ing to the geometric law with parameter ' 
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This result based on the memoryless property of the exponential distribution 
will be proved in Section [5711 

Let m be a nonnegative integer. The event [N = m + l} implies that the first sites to 
be infected are in chronological order (0, 0), (1, 0), . . . , (m + 1, 0) and finally (0, 1); 
see Figure |7] This provides the first moves of particles in the TASEP ^ obtained 
by the Rost's coupling. Precisely, Pq overtakes Hq, . . . , H„,+\, thus at time G(0, 1) 
particle Pi overtakes Hq. To sum up, on the event {N = m + l], ^(G(0, 1)) is equal 
to the configuration rf", introduced in (|7]l. 
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i- 



P^ P2 Hq Pi Hi 



H„, HmiA Pq H, 



Jmi-2 



Figure 6: On {iV = m + 1} and at time G(0, 1), the TASEP obtained by the Rost's 
coupling is equal to the configuration rf. The - and + pairs defined in Section [1]2] 
respectively consist of particles Hq and Pi and particles H^+i and Pq. 

Since G(0, 1) is a stopping time, the strong Markov property implies 

Lemma 8. Conditionally to {N - nt + I], the shifted process ^{- + G(0, 1)) is the 
TASEP with initial configuration rj'". 



4.3 Competition interfaces 

Let us recall that C{s) is the set of sites z € N^ whose geodesic goes by the source 
s, for s e {(0, 2), (1,1), (2, 0)). The aim of this section is to define the borderlines 
between the clusters C(2,0) and C(l, 1), and between C(l, 1) and C(0,2). 

The - competition interface is a sequence (^,7)">o defined inductively as fol- 
lows: ip~ = (0, 0), (f'^ = (0, 1) and for « > 1, 



fn+l 



<p- + ihO) if ^,; + (l,l)eC(0,2), 

if- + (0, 1) if ifi- + (1, 1) e C(l, 1) U C(2,0) . 



(15) 



In an equivalent way, ^^^^ chooses among the sites (p~ + (1,0) and (p~ + (0, 1) the 
first to be infected. Moreover, it is easy to draw the competition interface (i^,7)«>o 
from a realization of clusters C(2, 0), C(l, 1) and C(0, 2). Indeed, ip~ is the only site 
{x,y) € W- such that x + y = n,{x+ l,y) belongs to C(l, 1) U C(2, 0) and {x,y + 1) 
to C(0, 2). So, the directed path {(fi~)n>Q well describes the borderline between the 
clusters C(l,l) and C(0, 2). 

In the same spirit, the borderline between C(2,0) and C(l, 1) is described by the 
+ competition interface. This is a sequence {<p^)n>o defined inductively by ^^ = 
(0,0), (/)| -(1,0) and for «> 1, 



fn+l 



(^: + (l,0) if ^,t + (l,l)eC(l,l)UC(0,2), 
^: + (0,l) if ^,t + (l,l)eC(2,0). 
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When the competition interfaces (</',t)n>() and (9J^)„>o meet on a given site z.q (see 
the right hand side of Figure [T]) then they coincide beyond that site zo which is the 
larger (with respect to the L^-norm) element of C(l, 1): 

min{« > 1, (^,7 = ^n) = max{x + y, {x,y) € C(l, 1)} . 

In particular, there is coexistence if and only if the two competition interfaces never 
meet: 

\fn > 2, ^- + ^l . 

4.4 From competition interfaces to tagged pairs 

Let e € {+, -). Consider the competition interface (^^)n>o and its continuous-time 
counterpart, the interface process (^f defined by 



vf > 0, <f{i) = Y^ ^i\Gi.^„)fi(^i^,m 



«>0 

Set 

Vf > 0, {l\i), r(t)) ■= (p%t + G(0, 1)) . 

By construction of (^~)„>o, (p~{f) is (0,0) until G(0, 1) and 0"(G(0, 1)) is (0, 1). 
Besides, on the event {N = m + \], the point (f)^{G{0, 1)) is known too. Assume 
this event satisfied. On the one hand, sites (2, 0), . . . , (m + 1, 0) are infected before 
(1, 1), . . . , (m, 1) which yields (^+ ^ (2, 0), . . . , ^+^j = (m + 1, 0). On the other 
hand, at time G(0, 1), neither site (m + 2,0) nor site (m + 1, 1) are still infected 
which means ^,"^^2 ^^ '^°'- y^*^ determined. In conclusion, 0"*'(G(0, 1)) is equal to 
(m + 1, 0). To sum up, on the event {N = m+ I], 

(/"(O), 7"(0)) - (0, 1) and (/+(0), 7+(0)) - (m + 1, 0) . (16) 

See also Figure |7] 

Let ^ be the TASEP obtained by the Rost's coupling and assume the event {N = 
m+ 1} satisfied. Thanks to LemmaO we know that the shifted process ^(- + G(0, 1)) 
is the TASEP with initial configuration if". Recall that, in ^(f +G(0, 1)), the position 
of the oo particle of the e-pair is denoted by H^{t) (Definition |2l). Denote also 
by P^it) the position of the 1 particle of the e-pair. Of course, for any time t, 
P^{t) = H^{t) + 1. Moreover, at time t - (and always on{N -m + 1)), 

(//"(0),P~(0)) = (-1,0) and (//+(0),/'+(0)) = (m+ l,m + 2) . (17) 

The next result links competition interface (0'^(f + G(0, l))r>o to the e pair [oo 1]. 
Precisely, the coordinates {I^{t), J^{t)) are given by the labels of particles oo and 1 
constituting the s pair at time t. 

Lemma 9. The following identities hold on the event {N = m -\- \]. For any ? > 
and E € {+, -}, 

{H%t), P\t)) = {Hi.(,){t + G(0, 1)), Pjs(,^{t + G(0, 1))) , (18) 

14 



, 





























1 2 m m+l 





Figure 7: The infected region at time G(0, 1) conditionally to {A^ = m + 1), deli- 
mited by the black broken line. The two black squares represent (p^{G{0, 1)) and 
(p'*'{G{0, 1)). Since cp"^^^ chooses the earlier infected site among cpf^ + (1,0) and 
iffi + (0, 1), the interface (p^(t) is always in a corner formed by the black broken 
line. Combining with the algorithm given at the end of Section 14.11 this justifies 
heuristically why 0^(?) corresponds to a pair [oo 1] in the TASEP ^ obtained by the 
Rost's coupling. 



and 



Moreover, for any t >0, 



H%t) - r{t) - r{t) 



H\t) = H-{t) « {i\t),j\t)) = {r{t),r{t)) 

o ^+(f + G(0, 1)) - 0~(f + G(0, 1)) . 



(19) 



(20) 



Recall that Tcoi is the time at which the tagged pairs collide (Definition |2l). 
Assume N - m + \ and T^ox < °°- Then, just before time Tcoi, the two tagged 
pairs in the TASEP ^(- + G(0, 1)) are side by side and their labels satisfy I~{t) = 
I*(t) - 1 and J~{t) = 7+(0 + 1 (thanks to ([TS])). Thus, at time Tcoi, the configuration 
• • • [cx5 l][oo 1] • • • becomes • • • cx)[cx) 1]1 • • • and thenceforward the two interfaces 
coUide (thanks to (l20l)): 

0+(reoi + G(0, 1)) = 0-(rcoi + G(0, 1)) . 

Actually, Tcoi + G(0, 1) is the time at which the last site of C(l, 1) is infected. Fi- 
nally, remark the correspondence between competition interfaces and tagged pairs 
still holds after their collision. 
Lemma |9] will be proved in Section [ 



5 Proofs 

5.1 Proof of Lemma |7] 

Let m > 2 be an integer. First, 

P(A^ > m I A^ > 1) - P(A^ >m\N>m-l)x ¥{N > m - \ \ N > I) 
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By the memoryless property of the exponential law, 

P(A^ >m\N >m- I) 

= P(w(l, 0) + . . . + aj{m, 0) < co{0, 1) | w(l, 0) + ... + co(m- 1,0) < co(0, 1)) 
= P(w(m,0)< w(0, 1)) 

= 1/2. 

Hence, by induction we get P(A'^ > m \ N > I) = 2"™^' which is also true for 
m = I. This means that, conditionnally to {A'^ > 1), A'^ is geometrically distributed 
on {1, 2, . . .} with parameter 1/2. In other words, 

P(A^ ^ m I A^ > 1) = 2"'" (m > 1). 

5.2 Proof of Lemma |9] 

Throughout this proof, we assume N = m + \. Let us start with proving (fTSl) 
in the case e - -. In order to lighten formulas, let us denote by t„ the time 
G((^;;) - G(0, 1). Since ip^ = (0, 1), ti is equal to 0. At that time, 

{H'{0),P~m = (-1,0) 

= (//o(G(0,l)),Pi(G(0,l))) 

- (///-(0)(G(0,1)),P^-(0)(G(0,1))) , 

thanks to relations (IT6l) and ([TT] ). So, (ITSl) holds at time ti (and for e = -). Let us 
proceed by induction on times (t„)„>i. Assume (fTSl ) holds at time t„ for a given 
« > 1, i.e. I~{t„) and J~{Tn) are the labels of particles oo and 1 of the - pair at time 
T„, and prove it still holds for any time t € [r„; t„+i]. By definition, {I'{Tn), J~{Tn)) 
are the coordinates of the competition interface 0~(t„ +G(0, 1)) = (p~iG{(p~)) = ip~. 
At the next step, ^"^^ chooses the earlier infected site among (/~(t„) + 1, 7^(r„)) 
and {I^{Tfi), J^{Tn) + 1), say for example 

(r(T„+i), /-(r„+i)) = 0"(T„+i + G(0, 1)) = ip-^^ - {r{Tn) + 1, 7-(T„)) . 

Then, at time t„+i + G(0, 1) = G(i^~^j), particles Hi-(r„)+i and Pj-{t„) exchange 
their positions while Hj-^t,,) and Py-(r„)+i have not yet done (see the Rost's rule 
(fT4]) ). This statement has two consequences. The first move of the - pair after time 
T„ + G(0, 1) takes place at time t„+i +G(0, 1): (ITSl) holds for any time t € [t„; t„+i). 
Thus, at time t„+i + G(0, 1), the - pair jumps one unit to the right and its particles 
oo and 1 then become Hi-(^t„)+i and Pj-(r„)- So, 

(//-(T„+i), P-{Tn+,)) = (///-(r„)+l(T„+l + G(0, 1)), P/-(r„)(T„+l + G(0, 1))) 

= (^/-(r„,,)(^«+l + G{0, 1)), /'7-(r„,.)(^.+ l + G(0, 1))) , 

i.e. (fTSl) holds at time t„+i. The case (/j~^j = (/~(t„), J~{t„) + 1) leads to the same 
conclusion. 
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The case e = + is very similar. This time, put t„ = G{(p^^) - G{0, 1). We have 
akeady seen that at time G(0, 1) and on the event {N = m + I}, ^^^^^j = (m + 1,0) 
and ¥'m+2 ^^ ^^^ y^^ determined. So t„,+i < and 7^+2 > 0. Relation (fTSl) holds at 
time t = thanks to ([161) and (ITT] ): 

(//+(0), P+(0)) - (m + 1, m + 2) 

- (//„,+i(G(0,l)),Po(G(0,l))) 

- (///.(0)(G(0,l)),Py.(0)(G(0,l))) . 

Thus, the same induction as before, but on times (0, t^+z, Tm+3 • ■ •)> allows to con- 
clude. 

Let s € {+,-). It can be deduced from the previous remarks that when the e pair 
jumps one unit to the right, i.e. H'^ increases by 1, the label of its 00 particle in- 
creases by 1 whereas the one of its 1 particle remains the same. Conversely, when 
the e pair jumps one unit to the left, i.e. //^ decreases by 1, the label of its 1 particle 
increases by 1 whereas the one of its cx3 particle remains the same. To sum up, for 
any f, 

H%t) - H^'iO) = r{t) - r{t) - {r{0) - r{0)) . 

Combining with 

//-(O) = -1 - r(0) - 7'(0) and //+(0) = m+\= /+(0) - 7+(0) , 

we get ( fT9l ). 

It remains to prove (l20l ). Thanks to (fT9l) . the equality H^{t) = H'{t) is equivalent 

to 

r{t)-tit) = j-{t)-r(t). (21) 

Now, the directed character of the LPP model implies the differences /"(?) - /"""(f) 
and J^it) - /^(O are respectively nonpositive and nonnegative. So, (|2T1) forces 
I-(t) ^ l+{t) and J-it) = J+it). 

5.3 Proof of Theorem [1] 

In Section 1431 the coexistence phenomenon has been described in terms of com- 
petition interfaces: 



Vn > 2, «„ > o Wn>2, if^ i^ ip^ 



Let m be a nonnegative integer. Relation (EUb of Lemma |91 states that, on the event 
IN = m + \], the two competition interfaces {(p^)n>\ and ((/?^)„>i never meet if and 
only if the collision time Tcoi of the tagged pairs in the shifted process ^(- -1- G(0, 1)) 
obtained by the Rost's coupling, is infinite. Moreover, conditionally to {A^ = m-i- 1}, 
f(- + G(0, 1)) is the TASEP with initial configuration 77'" (Lemma[Hl). Let P^ be its 
probability measure. Then, 

P(Vn > 2, (^„ ^ ^+ I A^ = m + 1) - P^Crcoi - ~) • 
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The coupling stated in Section [331 between a TASEP with initial configuration 77'" 
and a three-type TASEP ^' with initial configuration ?7(3),m implies 

Pmcrcoi - ~) - p;(v?, 2[f (0] < 3[r(0]) , 

where P^ denotes the probability measure of ^'. Finally, the previous probability 
is equal to 1 - 2/(m + 3) (Proposition |4l). Combining the previous identities, it 
follows: 

P(V?i > 2, a„ > I A^ - m + 1) - 1 . (22) 

m + 3 

We conclude using symmetry of the LPP model, P(A'^ > 1) = 1/2, Lemma |7] and 



P(V?i > 2, a„ > 0) - 2P(V?i>2, a„ >0, A^> 1) 



= 2 2^ P(V?i >2, a„>0 , N ^m + \) 

m=Q 
00 

= 2 2 P(V« > 2, an>0\N ^m+\) 

m=0 



ZjOf^+A m + 3/ 

111— (I ' ' 



X P(A^ = m + 1 I A^ > 1) P(A^ > 1) 



m=0 

= 6-81og2 



The last equality comes from the formula 

log2^ V — 



, m2" 

m=l 

Let us point out here that, thanks to the memoryless property of the exponential 
distribution, initial conditions a>(0, 0) = w(l,0) = w(0, 1) = used in |[3l amounts 
to conditioning by the event {A'^ = 1}. So, their coexistence result (Theorem 4.1) 
corresponds to (l22l) with m = 0: 

P(V« > 2, a„ > I A^ - 1) - - . 

5.4 Proof of Theorem H 

Our goal is to prove that coexistence ahnost surely implies positive density: 

P(V«>2, Q'(«)>0 and lim — — = ) - . (23) 

y n->ca n I 

For E £ {+, -) and n > 1, let us denote by ^ the angle formed by the half-line 
[(0, 0), ipf^) with the axis y = 0: 
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Expressing the conditions a{n) > and lwaa(n)/n = in terms of angles 9^, 9^^ 
and using the symmetry of the LPP model with respect to the diagonal x - y,it is 
sufficient to prove 

P fv?i > 2, 9- > 9^, and lim 6*;; - 0+ = A^ > l) = 

\ ;l— >oo / 

or, in an equivalent way, that the conditional probability 

P fv« >2, 9-> 0,t and hm 6*; - 6*+ = A^ - m + l) (24) 

\ ;i— >oo / 

is null for any m e N. 

Let m be a nonnegative integer. In ||3, Ferrari and Pimentel have studied the 
asymptotic behavior of the border between the two subsets D(l, 0) and D(0, 1) of 
N^ formed by sites whose geodesic respectively goes by (1,0) and (0, 1). This 
border is described as a sequence {ipn)n>o -a competition interface- defined by 
(fo = (0, 0) and for n >0, 

^ r <p„ + {l,0) if ^„ + (l,l)eD(0,l), 
^"■•■^ \ ^„ + (0,l) if ^„ + (l,l)eD(l,0). 

When 6t»(l,0) < oj{0, 1) the geodesic of (1, 1) goes by (0, 1) rather than (1,0). In 
this case, 

D(l, 0) - {(1,0)} U C(2, 0) and D(0, 1) - {(0, 1)} U C(0, 2) U C(l, 1) . 

So the sequences (^„)„>o and (^^)n>o coincide on the event {N = m + 1} which is 
included in {// > 1} = {w(l,0) < w(0, 1)}. Now, by Proposition 4 of 15], (0+)„>o 
converges a.s. to a random angle 9. Thus, by Proposition 5 of IJH, 

hm = f{9) , (25) 

where / is a deterministic function (whose expression is without interest here). 
When the difference 9~ - 9'^ tends to 0, results of 1 5 ] apply again and yield (l25l) 
replacing + with -. Therefore, (l24l) is upperbounded by 

/ hm,^eo r' int) - r{t) - (/-(o - rm = 

\ and V?, /-(?) - /-(?) < /+(0 - r{t) "" 

Now, thanks to the Rost's coupling (Lemmas [8] and |9l relation (fT9l)). the above 
conditional probability is equal to 

/ H^{t)-H-{t) ^ \ 

Vm lim — — — = and V?, H'{t) < H^{t) , (26) 

where P„, denotes the probability measure of the TASEP with initial configuration 
77™. Finally, using Lemma[3j the quantity (l26l) becomes 

p;(iim 3[f(0]-2[r(0] ^ Q ^^^ w,2[r(0] < 3[r(0]l , 

\t^oo t ) 
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where £,' is a three-type TASEP with initial configuration ?y(^)"' and P'„ its proba- 
bihty measure. Lemma [5] achieves the proof of (l23l) . 

It remains to prove that a.s. the density of the cluster C(l, 1) cannot be equal to 
1. By symmetry with respect to the diagonal ;c = >', it suffices to show that 

P I hm — - 1 and w(l, 0) < oj(0, 1) | = . (27) 

When the density of the cluster C(l, 1) equals to 1, that of cluster C(2, 0) is null. 
In this case, the + competition interface (^,t)">o is asymptotically horizontal and 
the sequence (0^)„>o converges to 0. Furthermore, under the condition a»(l,0) < 
a<(0, 1), the competition interfaces (^,t)«>o and (^n)„>o -previously introduced in 
this proof- coincide. Proposition 4 of [5 1 then ensures the convergence almost sure 
of (0^)„>o to a random angle 6. To sum up, 

pf lim — = 1 and w(l,0) < w(0, 1)1 < P(6» = 0) . 

Theorem 1 of |5 ] also says the distribution of 9 has no atom. This proves (1271 ). 

It derives from the above arguments that cluster C(2, 0) has a positive density 
on the event {(x>(l,0) < a){0, 1)), i.e. with probability one half. Actually, this holds 
with probability 1 and the same is true for C(0, 2). To do so, let us remark that the 
cluster C(2, 0) grows when the weights 6t»(l, 0) and a>(0, 1) are exchanged, provided 
6ti(l, 0) is smaller than co(0, 1). It then can be proved that 

p(lim 0+ ^ and w(l,0) > w(0, 1)) < p(lim 6*+ ^ and aj{\,0) < aj{0, 1) 

\n—>cx3 J \n~>oo 

We have shown that the right hand side of the above inequality is null. Conse- 
quently, the probability of the event |lim 6'^j = 0) is null which implies that the 
cluster C(2, 0) has a.s. a positive density. 
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